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2" Abstract 

o\ • 

Properties of the neutron rich n Li nucleus are calculated in the framework of the 
cluster model 9 Li +n + n. The formalism of the harmonic oscillator representation is 
q ' used for the description of bound and continuum spectrum states in the three-body- 

democratic-decay approximation. It is shown that this approach allows one to take 
into account adequately the long asymptotic tail of the n Li wave function {neutron 
halo) and to reproduce correctly the binding energy, radius and n Li electromagnetic 
^ ■ dissociation cross-section on target nuclei. The shape and the energy position of 

the B(E1) peak corresponding to the soft dipole mode are also in agreement with 
experiment. 

Introduction 

Recently secondary beams of radioactive heavy ions become available. As a result one get 
a tool for experimental studies of light neutron-excess /3-unstable nuclei. n Li is one of the 
most interesting nuclei of the type, and its properties have being intensively studied both 



a 



U 



experimentally |I|-0 and theoretically (see, e.g., [f7|-[pl| and references therein). 

The main peculiarity of the n Li nucleus is the so-called neutron halo formed by two 
weakly-bound neutrons (two-neutron separation energy e 2 n = 247±80 keV; note, that there 
are no bound states in two-body subsystems ( 9 Li+n) and (n + n) p5|). As the binding 



energy of the neutron pair is small, the wave function decreases slowly at large distances. 
It results in the anomalously large value of the n Li r.m.s. radius, <r 2 > 11 = 3.16 ± 0.11 
fm (to be compared to the one of 9 Li, <r 2 >\ /2 = 2.32 ± 0.02 fm |25|). Evidently, the 
anomalously large electromagnetic dissociation (ED) cross section of n Li beam on heavy 
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target nuclei is the manifestation of the neutron halo (the n Li ED cross section on the Pb 
target is 20 times the one of the carbon beam Q). To explain this effect it was suggested 
[0, H that oscillations of the halo neutrons with respect to the core 9 Li give rise to the 
specific low energy branch of the giant dipole resonance (soft dipole mode) with excitation 
energy less then 1 MeV. Soft dipole mode is supposed to exhaust about 10% of the dipole 
energy- weighted sum rule (EWSR). 

The shell-model [0, P, [l|, H, cluster |]| 0, U, §§ II and combined cluster- 
shell-model [O, O, |24| approaches have been used in theoretical studies of the n Li neutron 



halo properties. In some of these calculations (see, e.g., (T(|-[|l||] a variational approach has 
been used which does not allow for continuum spectrum effects. There are also calculations 
with an allowance for two- 



24 continuum. 



PI |3 or three-body $Q [Lj| g|, [ 
In this contribution we present the study of the neutron halo properties of n Li in a 
three-body cluster model 9 Li+n + n. We make use of a three-body wave function expan- 
sion in six-dimensional harmonic oscillator eigenfunctions. Our attention is focused on the 
effects of three-body continuum. The continuum spectra effects are accounted for in the 
framework of the oscillator representation of the scattering theory p6j . This method has 
been generalized |27| on the case of the true three-body scattering |2^, |29| and has been 
successfully used recently in the study of the 12 C monopole excitations in the cluster model 
a + a + a [Q . As it has been noted, there are no bound two-body subsystems like ( 9 Li+n) 
or (n + n) in the system. Therefore, only the so-called "democratic" three-body decay 
[ f29| , |3J"H is inherent to this system. It can be adequately described in the framework of the 
"true" three-body scattering theory (democratic decay approximation) f29fl . In the study 



of the ground and excited states properties we go beyond the pure diagonalization of the 
Hamiltonian matrix on the basis of oscillator functions. We search for the S'-matrix pole 
corresponding to the ground state, and construct an infinite series expansion in oscillator 
functions for the ground and excited state wave functions. It enables us not only to calcu- 
late the n Li ground state energy with high accuracy, but to describe also the wide spatial 
distribution of two valence neutrons in n Li (neutron halo). 



1 The model 

The n Li ground state and the continuum spectrum wave functions are calculated in the 
framework of the three-body cluster model 9 Li+n + n. The model assumptions are the 
following ones: 

1) The cluster 9 Li is supposed to be structureless and the excitations of its internal 
degrees of freedom are not considered. 

2) We don't account for non-central components of the interaction between two valence 
neutrons and between valence neutron and the cluster 9 Li. Therefore, the wave function 
can be characterized by the three-body orbital angular momentum J and its projection M. 

3) The states with the total spin of the valence neutron pair S = are only considered, 
and the ground state three-body orbital angular momentum is supposed to be equal to 
zero: L = 0. 

4) n- 9 Li interaction is described by the shallow potential of Johansen et al |13| effec- 



tively simulating the Pauli principle j32|. iVTV-interaction is described by the Gaussian 



potential 13 



5) Only democratic decay channels are allowed for. 

The wave function of the system 9 Li+n + n, ^>jm(x, y), is expanded in three-body 
hyperspherical functions, & xl y JM (p), 

1>Mx,y) = E irf&i» ® l K lyJM (p) , (i) 

J\ L X Ly 

where K is hypermomentum, l x and l y are the angular momenta corresponding to the 
Jacobi coordinates 



^-(ri-r 2 ), y =y__(__-r 8 ), (2) 

respectively, m is the neutron mass, r, are coordinates of the valence neutrons (i = 1, 2) 
and the cluster 9 Li (i = 3), p= (x 2 + y 2 ) 1//2 is a three-body hyperradius, p stands for the 
set of angular variables in the six-dimensional space. 
In the c. m. frame the Hamiltonian is of the form: 

H = T + V 12 + V 13 + V 23 , (3) 

where T is the three-body relative motion kinetic energy operator, and V^- are the two-body 
potentials. For the radial wave functions ip K \ l (p) we have the usual set of the ii'-harmonic 
method coupled equations (see, e.g., |p§fl). The equations are solved by expanding the radial 
wave function, 



tfl l v (P) = E D n ] K l x l y (E)<PnK ( P ) , (4) 



n=0 

in the six-dimensional harmonic oscillator eigenf unctions 



Vuk(p) = (-l) n y r(n + ^ + 3) P K L« +2 (p 2 ) exp(-p 2 /2) , (5) 

where L"(x) is the associated Laguerre polynomial. 

For large values of the total number of oscillator quanta, N = 2n + K, the potential 
energy matrix elements are small compared to the matrix elements of kinetic energy which 
increase linearly with N. In the framework of the oscillator representation of scattering 
theory the matrix of the potential energy, V = V\ 2 + V13 + V23, on the basis of functions 
(Eh is truncated, i.e. all matrix elements V^ r for N > N and/or N' > N are neglected 
(JV denotes the truncation boundary). The kinetic energy matrix, T^ r r , is an infinite 
tridiagonal matrix for any value of the multi-index T = {K, l x , l y } which labels different 
decay channels in the hyperspherical representation. The eigenvectors of this matrix which 
are important for the description of the wave function asymptotical behavior at p — ► 00, 



in the oscillator representation are of the form ||27|| : 



D { n J) K fl( E ) = It ^rr' d!&{E) - 5 rr , C { J(E) ] for E> 0, (6) 



i(-0 fjp\ _„ r>(+) 



D":>(E)=a T C^{E) forE<0, (7) 



where Srr> is the matrix element of the S-matrix, and «r is the normalization constant of 
the bound state wave function for the channel I\ 

There are analytical expressions for the functions C^(E) 



i(±) 



C^(E) = C nK (E) ± iS nK {E) , (8) 



where 



(Hi 

and the total energy of the system, E = ^q 2 hu. The functions S u k(E) and C u k(E) give 
rise to functions with the following asymptotics in the coordinate space: 

oo -i 

£ S nK{E) <Puk(p) = / — - Jk+i(w) , (12) 

n=0 V^p 2 

oo 1 

£ C nJf (S) ^(p) ^ -7=^ A^ +2 (gp) , (13) 

n=0 Vft^p 2 

The asymptotic expressions (0) for Dl l j <lx ^ mma (E) correspond to the account of demo- 
cratic decay channels only [p9 , i.e. the wave function in the asymptotic region is a su- 



perposition of an outgoing six-dimensional spherical wave in the channel r" and ingoing 
six-dimensional spherical waves in all channels I\ The applicability of the democratic de- 
cay approximation for various nuclear reactions with few particles in a final state has been 



discussed in p9| , plf . In the framework of the oscillator representation of scattering theory 
this approximation has been used |3(| in the study of monopole excitations of 12 C in the 
cluster model a + a + a. The analysis of experimental data for 2 + states in A = 6 nuclei 
within the a + N + N cluster model assumptions has been performed using this approxi- 
mation in ref . |33|] . This analysis has demonstrated that the approximation is an adequate 
tool for description of the decay states of Barromean nuclei. 

The energy spectrum and the wave functions are calculated in the following way |27| . 
First of all, one should find eigenvalues, E ( x , and eigenvectors, {7„r }> °f the truncated 
Hamiltonian matrix, {H™ r }, 2n + K < N, 2n' + K' < N; i.e. one should solve an 
eigenproblem for the set of equations 

2n'+K'<N 

£ [ H nT {J) - <W ^rr' E^} 7^r' (J) = , 2n + K < N . (14) 

n', V 

In the usual variational approach, the lowest eigenvalue, Eq , and the corresponding 
eigenvector, J7„r ~ } are treated as the ground state energy and the corresponding wave 
function in the oscillator representation. Such approximation is unable to describe the 
slowly dying asymptotic tail of the halo neutron space distribution. Thus, in order to 



describe the neutron halo properties of the n Li nucleus, one is pushed to account for the 
asymptotic region N > N. For bound states this account is equivalent to the location of 
S-matrix poles. The S-matrix can be calculated by the equation [27] 

s = (# ) r i # ) . (15) 

The matrices A^ and A^ have the following matrix elements: 

4? = <Vr< Cgr(£) + Pnf(E) l£/ fl , +1 C% >K ,{E) , (16) 



the kinetic energy matrix elements, T% n+1 = —\huj J(n + l)(n + K + 3), n = -|(iV — K) 



n, n+1 
(AT — K'\ Q r,r1 



n' = UN -K'), and 



(A) (A) 

a E x - h 

To calculate the bound state energy, i.e. to locate the corresponding S'-matrix pole, one 
should solve the nonlinear equation 



det A {+) = . 



which can be easily obtained from (15). Asymptotic normalization constants of the bound 



states, ar, can be found by numerical solution of the set of linear homogeneous equations 

J2 4p «p = • (19) 

r' 

The set of «r can be multiplied by any common multiplier. This uncertainty is eliminated 
by numerical normalization of the wave function. Now the coefficients D nKl t (E) of the 

expansion @) in the asymptotic region N > N can be easily obtained using (|7|). 

For the continuum spectrum states we calculate S'-matrix for any positive energy E 
using (|T5l)-(|T7l), and than calculate the set of E] x ' Klxl (E) in the asymptotic region N > N 
using (0). 

The coefficients D nKlxl (E) in the inner region N < N are calculated by the expression: 

D nKi x i y ( E ) = _ zJ P wv(E) T Wfl , +1 D h , +ljKll ,^{E) , (20) 

r' 

where h' = |(iV — K') and for calculation of -D-+! Kix i (E) eqs. (|7J) or (|7J) are used for 
continuum or bound state, respectively. 



2 The results of the calculations 

The interactions of the valence neutrons with each other and with the cluster 9 Li are 
described by the potentials 1^12(^12) and Vi3(ri 3 ) = 1^3(^3), respectively. We use the 



following parametrization of the potentials [p 





Vij(r) = 


= V$> exp [ 


-(r/bff) 2 } 


+ \f expl-(r/bffy] 




T/(l)_ 


-31 MeV, 


V V2 


= 0, bty = 1.8 fm 


r(l) _ 
13 _ 


-7 MeV, 


V {2) - 


-1 MeV, 


&&> = 2.4 fm, fcg> 



3.0 fm. 

In the external asymptotic region N > N we consequently allow for channels T char- 
acterized by K = K min , K min + 2, . . . (K min is the minimal possible value of K for a 
given J) until the convergence for all physical properties under consideration is achieved. 
The convergence is found to be very good, and the allowance for the decay channels with 
K > K min + 2 do not yield any visual variation of the results. So, we consider in the external 
asymptotic region N > N the channels with K < K min + 2 only. Note, that components 
with all possible values of K < N are accounted for in the calculation of the wave function 
in the inner region N < N. 

The parameter Twj is set to be equal to 7.1 MeV in our calculations. This value corre- 
sponds approximately to the minimum of ground state energy E . 



2.1 The ground state 

The results for the n Li ground state for different values of the truncation parameter N 
are presented in the table 1. The variational ground state energies, Eq , obtained by the 
pure diagonalization of the truncated Hamiltonian matrix are listed in the second column, 
while the J-matrix results, E , which are the solutions of the eq. fll~8|) are listed in the third 
column. It is seen, that by locating the S'-matrix pole using eq. (|i~8|) that is equivalent to 
the allowance for the long asymptotic tail of the wave function, we improve essentially the 
convergence for the binding energy. 

The results presented in the table 1 have been obtained using Lanczos smoothing of the 
three-body potential energy matrix [34]. It has been shown in ref. [34] that the smoothing 



improves the convergency of calculations of two-body scattering phases. The results of the 
full J-matrix calculations of the ground state energy of n Li obtained with the use of the 
smoothing, Eq, and without it, Eq ' s , for various values of the truncation boundary N are 
presented in table 2. It is seen from the table 2 that the smoothing causes underestimation 
of the binding energy of a three-body system. Nevertheless, Eq' s shows a staggering 
as N increases^. The smoothing results in a usual-type monotone decrease of E as iV 
increases and in a better convergence of continuum spectrum calculations. At the same 
time, the ground state wave functions obtained with the smoothing and without it differ 
very slightly. So, below we shall discuss only the results of calculations with the smoothed 
potential energy matrix. The role of the smoothing in three-body calculations we shall 



discuss in more detail elsewhere 22 



The n Li r.m.s. radius, < r 2 > 1 { , can be calculated by the following equation: 

< r 2 > n = -— < r 2 > 9 +- < p 2 > , (21) 

11 lima; 



'■Note, that J-matrix calculation is not a variational one, so, the ground state energy may behave 
non-monotonically as N increases. 



Table 1: n Li ground state properties (see text for details). 



Truncation 
boundary iV 


Ground state energy, 
MeV 


Neutron halo 

mean square radius 

< r 2 >}{ 2 , fm 


E {d) 


E 


< r 2 >if W 


< r 2 >l( 2 


12 
16 
20 

24 


-0.012 
-0.116 
-0.171 
-0.202 


-0.150 
-0.199 
-0.225 
-0.240 


2.83 
2.91 
2.98 
3.04 


3.31 
3.29 
3.31 
3.32 


Experiment 


-0.247±0.080 


3.16±0.11 



Table 2: The effect of the Lanczos smoothing on the ground state energy of n Li (see text 
for details). 



N 


10 


12 


14 


16 


18 


20 


22 


24 


E 


-0.119 


-0.150 


-0.181 


-0.199 


-0.215 


-0.225 




-0.240 


J^(W.S.) 


-0.136 


-0.266 


-0.204 


-0.262 


-0.232 


-0.262 


-0.249 


-0.263 



where < r 2 > 9 is the 9 Li r.m.s. radius and the mean square value of the hyperradius, 
< p 2 >, can be easily calculated using the ground state wave function. The values of < 
r 2 > 1 / 2 w and < r 2 > x l 2 obtained by the pure diagonalization of the truncated Hamiltonian 
matrix and with the allowance for the asymptotic tail of the wave function, respectively, 
are presented in the 4-th and the 5-th columns of the table 1. 

Figure 1 presents the transverse momentum distribution of the cluster 9 Li in n Li. This 
momentum distribution is currently supposed (see, e.g., |16|) to be proportional to the 9 Li 
transverse momentum distribution J^- in the fragmentation of high-energy n Li beams on 
target nuclei. The experimental data for n Li fragmentation [0 are also presented on the 
figure. The results of J-matrix calculations with N — 16 and N = 24 are so close to each 
other that their plots on fig. 1 coalesce and cannot be distinguished. 

It is seen that in calculation of the ground state, the allowance for the wave function 
asymptotics is very important for a weekly-bound system like 11 Li. The terms of expansion 
(fD with the number of total oscillator quanta N ~ 100 that cannot be obtained in the usual 
oscillator-basis variational calculations, play an essential role in the formation of the trans- 
verse momentum distribution, r.m.s. radius, etc. The convergence of < r 2 > 1 / 2 , transverse 
momentum distribution and other properties of the wave function (e.g., of the weights of 
its components) in the full J-matrix calculation is very good. Nevertheless, it is seen that 
the r.m.s. radius converges to a value that is somewhat larger than the experimental one, 
and the calculated transverse momentum distribution is narrower than the experimental 
one. These shortcomings can be overcome by the adjustment of n- 9 Li potential. We have 
not aimed to fit the potential to the n Li properties, we have just take its parameters from 
ref. ]n|. The differences between the calculated and experimental values of r.m.s. radius 
arise from the fact that in ref. |13| the potential has been fitted to the n Li r.m.s. radius 



and binding energy in variational calculations that are unable to reproduce these quantities 
with high accuracy. 



Figure 1: The 9 Li transverse momentum distribution in the ground state of n Li. 1 — vari- 
ational calculations, 2 -- J-matrix calculations. Experimental data for the 9 Li transverse 
momentum distribution in the fragmentation of n Li 800 MeV/nucleon beam are taken from 
ref. 01. 



2.2 The soft dipole mode 

The dipole transition operator in our model is of the form 

N 7 
M(Eln) = --^eyY llt (H), (22) 

where e is the proton charge, A = 11, Z = 3 and the number of valence neutrons, N v = 2. 
The operator (B^) corresponds to the excitation of the three-body cluster modes only. The 
excitation energy of the first excited state of 9 Li is relatively high (~4 MeV). So, low- 
energy .El-transitions correspond to the excitation of the cluster degrees of freedom only 



and should be described by the operator fl2"2"|) . 



The cluster reduced probability of the El-transition, B(E1; Ej — Eq), associated with 
the operator (p2|), is displayed on the figure 2. In J-matrix approach, zero- width pure 
variational peaks of B(E1; Ef — E ) gain finite widths and shift to lower energies. It is seen 
that in calculation of B(E1\ Ef — Eq) it is important to allow for the asymptotics not only 
for continuum states, but for the ground state as well. Due to the wide spatial distribution 
of the halo neutrons density, large distances contribute essentially to the El strength in the 
low-energy region. To account for this contribution and to get a convergence, one should 
allow for ground and continuum state wave function components with very large values of 



Figure 2: Cluster B(El;g.s. — > continuum) for n Li. Zero-width peaks with diamonds on 
the top (1) given in arbitrary units have been obtained in the pure variational approach 
for both ground and exited states, i.e. the wave functions for the initial and the excited 
states have been obtained by the diagonalization of the truncated Hamiltonian matrix. The 
curve with dots (2) has been calculated with the variational ground state wave function, 
the continuum spectrum wave function has been obtained by the J-matrix approach. The 
solid curve (3) is the result of J-matrix calculations for both the ground and continuum 
spectrum states. 

the total number of oscillator quanta N ~ 2000| in calculations of B(E1; Ef — Eq) . The 
large-distance .El-transitions enhance the peak of the cluster B(E1; Ef — Eq) and shift it 
to a lower energy. Obviously, this peak should be associated with the soft dipole mode. 
Thus, the neutron halo manifests itself in the appearance of the soft dipole mode that arises 
from the shift and enhancement of the B(E1; Ef — E Q ) peak. This effect appears to be 
very important in the calculation of the n Li electromagnetic dissociation cross section (see 
below) . 

Large distances contribute essentially to E2- and EO-transitions, too. The shift and 
the enhancement of B(E0; Ef — Eq) and of B(E2; Ef — Eq) are even more pronounced. 
Nevertheless, El- and EO-transitions do not play an important role in the electromagnetic 
dissociation, and we shall not discuss them here. 

Figure 3 shows the comparison of the results of our calculations of cluster B{E\\ Ef—E ) 
with the parametrization of experimental data of ref . . The agreement is reasonable. The 
form of the B(E1; Ef — E ) peak is well reproduced, the discrepancy in the position of the 
B(E1; Ef — E Q ) maximum is supposed to be eliminated by the adjustment of the potentials. 
The results of the B(E1; Ef — E Q ) calculations of refs. |L9|, |23| are also depicted on fig. 3. 
The three-body cluster calculations with the allowance for democratic decay channels of ref. 

§The classical turning point for the basis function with the total number of oscillator quanta N — 1000 
is at a distance of w!08 fm from the origin. 



Figure 3: Comparison of our results for B(El;g.s. — > continuum) in n Li with results 
of other authors. 1 - - this work (J-matrix method), 2 - - ref. |2"5[] , 3 - - ref. [I9~|, 4 - 
experimental data parametrization of ref. ||. 



nicely reproduce the energy of the soft dipole mode. Nevertheless, the form of the peak 
in our calculations is reproduced better. Note, that the authors of ref. p3[ used another set 
of the potential parameters. The calculations of ref. [||J with the allowance for two-body 
decay channels failed to reproduce both the position and the form of the B(E1; Ef — E ) 
peak. 

The soft dipole mode exhausts about 90% of the cluster EWSR associated with the 
operator 
that 



2j), iS^™.*. (£71). Nevertheless, using the results of ref. |35| it is easy to obtain 



>S clust( E1 ) _ * ~cq% 



(23) 



where S^ Q ^(E1) is the total EWSR accounting for excitations of all nucleons. So, the 
contribution from the soft dipole mode to the total EWSR is relatively small. In the 
vicinity of the sharp B(E1; Ef — E ) maximum at the excitation energy E «l-2 MeV 
only ~8% fraction of the total EWSR is exhausted. Nevertheless, the account for the soft 
dipole mode results in an essential increase of the electromagnetic dissociation cross section 
of 0.8 GeV/nucleon n Li beams on Pb and Cu targets. The wide space distribution of 
the halo neutrons density is also well-manifested in the electromagnetic dissociation of n Li 
beam. The contribution of the large-distance £7A-transitions to the cross section is about 
50%. The calculations have been performed by the equivalent photon method 
only parameter of the method is the minimal value of the impact parameter b 
for b m in the values of 9.0 fm for Pb and 6.8 fm for Cu target nuclei, respectively. These 



36 



The 

„„„. We use 



quantities are the sums of the n Li and target nucleus charge radii. With these values of 
b m in we obtain for the electromagnetic dissociation cross sections the values of 0.966 barn 
for the Pb target and 0.132 barn for the Cu target; the corresponding experimental values 
are 0.890 ± 0.110 barn and 0.21±0.04 barn, respectively [|J. E0- and ^-transitions give 
only 1.2% contribution in the cross sections. 

El-transitions in n Li have been studied in the framework of RPA+two-body-continuum 



model in refs. |15| ^(J. Though our model assumptions differ significantly from the ones of 
refs. [|15|, [2(J], the results are in good agreement. For example, the excitation energy values 
corresponding to the peaks of the function B(E1; Ef — E ) displayed on the fig. 2, are 



very close to the values that one can find in refs. |TR pQ 



Conclusions 

It is shown, that cluster model 9 Li+n+n yields a good description of the ground state prop- 
erties and El-transitions in the n Li nucleus. The oscillator functions expansion technique 
may be used in the studies of weakly-bound systems with long-tailed wave functions, e.g., 
in the study of neutron halo properties. For both bound and continuum states the correct 
account of the wave function asymptotics in the framework of the oscillator representation 
of scattering theory is very important in such studies. Low-energy El-transitions in n Li 
are of the cluster nature. The widths and the position of resonant states calculated in the 
democratic decay approximation are in a reasonable agreement with experiment. 

Acknowledgments 

We are thankful to Profs. J. Bang, B.Danilin, I.Thompson and J.Vaagen for valuable dis- 
cussions. 



References 

[1] Tanihata Let al. Phys. Rev. Lett. 1985, 55, 2676. 

[2] Kobayashi T. et al. Phys. Rev. Lett. 1988, 60, 2599. 

[3] Tanihata I. et al. Phys. Lett. 1988, B206, 592. 

[4] Kobayashi T. et al. Phys. Lett. 1989, B232, 51. 

[5] Anne R. et al. Phys. Lett. 1990, B250, 19; Riisager K., Anne R., Arnel S.E. Nucl. 
Phys. 1992, A540, 365. 

[6] Sackett D. et al. Preprint MSUCL, Michigan State University, 1993. 

[7] Ikeda K. INS report JHP-7, 1988 (in Japanese). 

[8] Tanihata I. Nucl. Phys. 1990, A520, 411c. 

[9] Bertsch G.F., Foxwell J. Phys. Rev. 1990, C41, 1300. 



[10] Hayes A.C., Strottman D. Phys. Rev. 1990, C42, 2248; Hayes A.C. Phys. Lett. 1991, 
B254, 15. 



[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 

[20 

[21 
[22 
[23 

[24' 
[25 
[26 

[27; 

[28; 
[29 

[30 

[31 



Hashimo O., Sagawa H., Arima A. Nucl. Phys. 1991, A523, 228. 

Tosaka Y. et al. Progr. Theor. Phys. 1990, 83, 1140. 

Johansen L. et al. Phys. Lett. 1990, B244, 357. 

Suzuki Y., Tosaka Y. Nucl. Phys. 1990, A517, 599. 

Fayans S.A. Phys. Lett. 1991, B267, 443. 

Zhukov M.V. et al. Nucl. Phys. 1991, A529, 53. 

Danilin B.V. et al. Preprint CNP-92/17, University of Surrey, 1992. 

Bang J.M., Thompson I.J. Phys. Lett. 1992, B279, 201 

Bertsch G.F., Esbensen H. Ann of Phys. 1991, 209, 327; Esbensen H., Bertsch G.F. 
Nucl. Phys. 1992, A542, 310. 

Lenske H. Proc. Int. Conf. Nucl. Struct, and Nucl. React, at Low and Intermed. 
Energ., Dubna, 1992, p. 132. 

Thompson I. J., Zhukov M.V. Preprint CNP-93/9, University of Surrey, 1993. 

Lurie Yu.A., Shirokov A.M., Smirnov Yu.F. To be published. 

Bang J.M., Danilin B.V., Thompson I.J., Vaagen J.S., Zhukov M.V. Private commu- 
nication. 

Zhukov M.V. et al. Phys. Rep. 1993, 231, 151. 

Ajzenberg-Selove F. Nucl. Phys. 1988, A490, 1. 

Yamani H.A., Fishman L. J. Math. Phys. 1975, 16, 410; Filippov G.F. Yad. Phys. 
1981, 33, 928; Nechaev Yu.L, Smirnov Yu.F. Yad. Phys. 1982, 35, 1385. 

Smirnov Yu.F., Shirokov A.M. Preprint ITP-88-47P, Institute of Theoretical Physics, 
Kiev, 1988. 

Merkuriev S.P., Faddeev L.D. Quantum scattering theory for systems of few bodies. 
Moscow, Nauka Publishes, 1985 (in Russian). 

Jibuti R.I., Krupennikova N.B. Hyperspherical harmonics method in quantum me- 
chanics of few bodies. Tbilisi, Metsniereba, 1984 (in Russian). 

Mikhelashvili TYa., Shirokov A.M., Smirnov Yu.F. J. Phys. 1990, G16, 1241. 

Jibuti R.I Elem. Part, and Atom. Nucl. 1983, 14, 741. 



Jensen A.S., Riisager K. Nucl. Phys. 1992, A537, 45. 

Bochkarev O.V. et al. Nucl. Phys. 1989, A505, 215. 

[34] Revai J., Sotona M., Zofka J. J. Phys., 1985, Gil, 745; Mares J. Czech. J. Phys., 
1987, B37, 665. 

[35] Sagawa H. et al. Nucl. Phys., 1992, A543, 575. 

[36] Winter A., Alder K. Nucl. Phys. 1979, A319, 518. 



